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Abstract. The energy of a Gamma-Ray Burst is one of the most in- 
teresting factors that can help determining the origin of these mysteri- 
ous explosions. After the discovery that GRBs are cosmological it was 
thought, for a while that they are standard candles releasing ~ lO^^ergs. 
Redshift measurements that followed the discovery of GRB afterglow re- 
vealed that GRBs are (i) further than what was initially believed and (ii) 
have a very wide luminosity function. Some bursts revealed an enormous 
energy release with a record estimate of 1.4 x lO^^ergs for GRB990123. 
The energy budget was stretched even further when it was suggested that 
the conversion efficiency of producing gamma-rays is low. The realization 
that GRBs are beamed changed this perspective and brought the energy 
budget of GRBs down back to a "modest" ~ 5 x 10^'^ergs. I discuss here 
various estimates showing that GRB energy is narrowly distributed and 
discuss the implications of this conclusion to GRB models. 



1. Introduction 

The energy release in a Gamma ray burst is one of the best clues on the nature 
of these objects. However, this simple and basic quantity is rather hard to 
get. The observed flux was readily available even for the first bursts. But 
the lack of a reliable distance estimate was a series obstacle. When redshift 
measurments became available it turned out that even this knowledge was not 
sufficient. Other complications arose. On one hand theoretical considerations 
showed that the efficiency of conversion of the energy to 7-rays could not be 
very high and estimates based on 7-rays alone were too low. On the other 
hand, it turned out that GRBs are beamed and estimates that assumed isotropic 
emission were too high. In this lecture I discuss the current understanding of 
the energetics of GRBs and the implication to models of the "inner engines" 
that power GRBs. 

In the first GRB revolution BATSE (Meegan et al., 1992) demonstrated 
that GRBs are cosmological. This has changed the distance scale by six orders 
of magnitude and the energy scale by twelve orders of magnitude. BATSE's 
count distribution for long bursts (Tgo > 2sec) is consistent with a cosmological 
standard candle distribution (Cohen & Piran, 1995) with E ~ lO^^ergs. This 
has led to the believe that (at least long) GRBs are standard candles and that 
the bursts are observed by BATSE up to z 2. 

The second GRB revolution took place in 1997 with the discovery of GRB 
afterglow (Costa et al, 1997; van Paradijs et al., 1997). The Italian-Dutch 
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satellite BeppoSAX provided angular position of several dozen long bursts to 
within about 3 arc-minutes which enabled follow up observations in the x-ray 
(see Piro, 2000), optical, milli- meter and radio frequencies which has provided 
a wealth of information on these explosions. In several cases the redshift of the 
afterglow or the host galaxy could be measured. This provided a new and direct 
estimate of the distance and of the energy involved. The first redshift, z = 0.835 
was obtained for GRB970508. The corresponding energy, 5.4 x lO^^ergs, was 
more or less in line with the standard candles estimates! 

It quickly became apparent that GRBs are not standard candles. With a 
redshift 3.4 and energy 2 x lO^^ergs GRB971214 was too far and too energetic 
for the standard candle picture. GRB990123 was even more energetic. It turned 
out that the GRB energy distribution is very broad and that their rate follows 
the star formation rate ( Totani, 1997; Sahu et al, 1997; Wijers et al, 1998). The 
consistency of a standard candle cosmological model with BATSE's peak counts 
distribution remains an inexplicable coincidence. 

The observations of extremely large energies (at least in some bursts) sug- 
gested that some GRBs involve energy release of stellar mass or more, ruling out 
several of the leading models at the time. The situation was even more worri- 
some when one considered the issue of efficiency. Most models and in particular 
the internal shocks model, cannot produce 7-rays at 100% efficiency. There was 
a concern that the efficiency of the internal shocks model is rather low (at the 
level of a few percent (Kobayashi et al, 1997; Daigne, &; Mochkovitch, 1998). 
This increased further the energy requirement for GRBs, ruling out most models 
and leading to a GRB energy crisis (Kumar, 1999). 

On the other hand GRBs could be beamed. The energy budget would then 
drop down significantly, by a factor 0^/2, where 6 is the opening angle of the 
jet. Rhoads (1999) pointed out that an expanding relativistic jet would exhibit 
spherical like evolution as long as F > 9~^. It will expand rapidly sideways 
afterwards. This behaviour will produce a break in the afterglow light curve. 
Using this break one could determine the opening angles and estimate the true 
energy budget. Already in 1999 Sari, Piran &; Halpern pointed out that the two 
most energetic GRBs (at that time) exhibit jet-like behaviour. This suggested 
that beaming is a strong factor in the energetics of GRB. 

I summarize here several different estimates to the energy of GRBs. I argue 
that there are good indications that different GRBs emit relatively constant 
energy. GRBs are standard candles after all. However, differences in beaming 
factors lead to a variation of three orders of magnitude in the isotropic equivalent 
energies and to a very wide apparent luminosity function. 

2. The Overall Picture and the Fireball Model 

Following the discovery of GRB afterglow we have learned during the last four 
years a great deal about long duration gamma-ray bursts (GRBs). These ob- 
servations are described well by the relativistic fireball model (sec e.g. Piran, 
1999). According to this model the energy from the central source is deposited 
in material that moves very close to the speed of light. The kinetic energy of 
this material is converted to the observed 7-rays as a result of collisions between 
fast moving material that catches up with slower moving ejecta. The afterglow 
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Figure 1. The internal-external shocks fireball model: The GRB is produced 
by internal shocks within the flow. The afterglow is produced by shock heated 
circum-burst matter. 



is produced latter by the shock heated circum-burst medium (see Fig. 1). The 
nature of the "inner engines" that expels relativistic material, which is respon- 
sible for the the GRBs, is not determined yet. I discuss here the energetic of 
GRB sources in the context of this model. One of the problems in discussing 
the energy of GRBs is that there are several different energies involved. I begin 
with a discussion of the different energy definitions.. 

2.1. Eneregy - Theoretical Considerations 

Our goal is to find the best estimate for 

• Etot- the total energy emitted by the source. 

However, it is practically impossible to determine this energy. It is not known 
what are all forms of energy release by the source. For example, the fireball 
model is based the conversion of the kinetic energy of relativistic flow (particles 
or Poynting flux) to 7-rays. It is impossible to detect non-relativistic particles 
that may also be emitted by the source. A second, more accessible quantity is 

• Erei- the energy of the relativistic flow emitted by the source. 

A fraction, e, of Erei is converted via the internal shocks to the observed prompt 
7-ray emission. The rest, (1 — e)Ej.ei, is dissipated later via external shocks on 
the circum-burst matter producing the afterglow. An unknown fraction of this 
energy is dissipated in the radiative phase (which has not been observed yet) 
during the first half hour of the afterglow. The remaining energy is: 

• Ek- the kinetic energy during the adiabatic afterglow phase. 

Ek is dissipated gradually over a period of months or even years producing the 
observed afterglow. Clearly, we have: Ek < (1 — e)Erei < Etot ■ Observations 
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of long time tails of GRBs suggest that this losses arc small (Burenin et al, 
1999, Gibhn et al, 1999; Tkachenko et al., 2000). Thus, Ek {I - e)E^^i. The 
efficiency, e, cannot be too large otherwise there won't be afterglow. It cannot 
be too small either, otherwise we will reach a GRB energy crisis. Combining 
these facts we expect, therefore, that Ek ~ Eyei to within say a factor of a few. 

2.2. Eneregy - Observational Considerations 

Given the observed 7-ray fiuence and the redshift one can easily estimate 

• E^^iso- the energy emitted in 7-rays assuming that the emission is isotropic. 

E^^iso can also be estimated from the BATSE catalogue by fitting the flux dis- 
tribution to theoretical models (Cohen &: Piran, 1995; Schmidt, 2001). As af- 
terglow observations proceeded, alarmingly large values (Kulkarni et al. 1999) 
(3.4 X lO^^'ergs for GRB990123) were measured for E^^iso- As wc discuss shortly 
(Rhoads, 1999; Sari Piran &; Halpern, 1999) GRBs are beamed and E^^iso is far 
larger than the actual energy emitted in 7-rays. We define instead: 

• E-y = {(P' l'2)E^^iso- 

Here Q is the effective angle of 7-ray emission. 

One would expect that E^ is a good estimate to E^^i as: £^-y = eErei- 
However there are several prooblems. First e is unknown. Moreover one would 
expect it to vary significantly from one burst to another. Second, the large 
Lorentz factor during the 7-ray emission phase, makes the observed E^ rather 
sensitive to angular inhomogeneities of the relativistic ejecta (Kumar Sz Piran, 
2000). During the GRB phase the relativistic Lorentz factor is at least a few 
hundred (See e.g. Lithwick,&; Sari, 2001 for a summary of the arguments con- 
cerning pair production opacity). Thus, the observed 7-rays come from a region 
whose angular size is 7"^ < 10~^. This is narrower by a factor of ten than the 
angular width of the narrowest observed jets. Thus the observed 7-rays span 
only a small fraction of the actual emitting region. The estimated energy based 
on this data alone could be misleading. 

2.3. The EflRciency of Internal Shocks 

The conversion efficiency of kinetic energy to 7-rays, e, depends on two factors: 
the conversion efficiency of bulk kinetic energy to energy of accelerated electrons 
and the efficiency of radiating this energy to 7-rays. It is usually assumed that 
this second factor is close to unity. The first factor depends on the strength of the 
relevant shocks and in turn on the relative Lorentz factors between the different 
shells (Kobayashi, Piran & Sari, 1997). There have been numerous attempts to 
estimate this efficiency (Kobayashi, Piran &: Sari, 1997; Daigne, &: Mochkovitch, 
1998; Kumar, 99; Spadaet al, 2000; Beloborodov, 2000; Kobayashi, & Sari, 2001; 
Guctta et al, 2001) and the results range from a few percent to near unity. The 
conversion can be efficient (close to unity) if the distribution of Lorentz factors is 
very wide (Beloborodov, 2000, Kobayashi, Sz Sari, 2001). In those cases in which 
GRB afterglow is observed this efficiency could not be too large. Otherwise there 
would be no energy left to produce the afterglow. 
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Figure 2. From Granot et al, (2001) Left: Density countor and velocity fields 
in an expanding jet long after the jet break. Right: Light curves from a full 
hydrodynamic calculation of an expanding jet. A clear break is seen, approximately 
at the time given by Eq. 1. 

2.4. Jets and Beaming 

The original fireball model assumed spherical symmetry. As always this was 
partially because of simplicity. However, the spherical approximation is valid as 
long as r > A9~^, where A9 is the scale of angular inhomegeneity. Because of 
time dilation sideways propagation is too slow and the information on the inho- 
mogeneity could not propagate (Piran, 1995). As the blast wave is slowed down 
by the circum burst matter the Lorentz factor decreases and when F ~ theta~^ 
it begins to expand sideways (Rhoads, 1999). This produces a break at the light 
curve, which is accompanied by a change in the spectral index (Sari Piran &; 
Halpern, 1999). Various models (usually analytic or semi-analytic) have been 
used to describe the hydrodynamic evolution of this expanding phase (Rhoads 
1999; Sari Piran &: Halpern, 1999; Panaitescu & Meszaros 1999; Moderski, Sikora 
& Bulik 2000; Kumar & Panaitescu 2000). Different assumptions have lead to 
different relations between the opening angle and the tb, the time of the break 
in the light curve. Following Sari Piran & Halpern (1999) I adopt here : 



where Ek,5i is the adiabatic kinetic energy in units of 10^^ ergs, -E51 = 2Ek,5i/0'^ 
is the isotropic-equivalent kinetic energy and tb^days is the break time in the 
afterglow light curve expressed in days. 

Detailed numerical simulation (Granot et al, 2001) show that the sideway 
expansion is more complicated than what was assumed in the simple analytic 
models (Fig. 2 depicts the density and the velocity field from a jet long after 
the jet break. The sideway expansion is not as prominent as expected.) Still, 
somewhat surprisingly, a jet break arises more or less at tb according to Eq.l. 

3. Energy Estimats 

The simplest and most direct estimates are of -E-y^jso- Schmidt (1999) have 
fitted the {V/V)max) distribution of BATSE's GRBs a cosmological model that 
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follows the star formation rate. He finds a break energy of 1.2 x 10^^ ergs in 
the 10-1000 keV band but as the higher slop is steeper the average energy is 
2.4 X lO^^ergs. This is of course the isotropic energy. The energy estimate 
depends on the specific cosmological model (Schmidt mentions variation by a 
factor of 2.5 between different models), on the star formation rate and on the 
assumptions on the average spectral indices of the GRB. 

For GRBs with known rcdshifts one can estimate E^^iso- Here care should 
be taken to obtain an exact estimate of the spectrum (sec Jimenez et al., 2001) 
and to add a proper K correction (Bloom et al., 2001). Fig 3 depicts the isotropic 
energy and the redshift for 17 bursts with known afterglow from Bloom et al., 
(2001). One sees a large spread (as seen also in Schmidt's (1999) luminosity 
function). Jimmenz et al (2001) consider 8 BATSE bursts with a spectroscopic 
redshift estimate and four additional bursts whose redshift has been estimated 
on the basis of the host galaxy R magnitude. They find an average E^^iso of 
1.3 X lO^^ergs. This agrees with the average of Bloom's sample but is higher by 
a factor of 5 than Schmidt's average for the full BATSE sample. 

The next step would be to correct E^/i^o for beaming. Frail et al. (2001, 
hereafter Fl) estimated E^ for 18 bursts with redshift. They find a very narrow 
distribution with typical values around 5 x 10^°ergs (see Fig 3) and FWHM of 
a factor of 5. This energy estimate should be taken with care as it depends crit- 
ically on the estimated jet opening angles di. Frail et al (2001, hereafter FOl) 
estimate the opening angles using Eq. 1. Panaitescu & Kumar (2001, hearafter 
PKOl) perform a detailed modeling of the afterglow and obtain different esti- 
mates for 9i. A comparison of the opening angles obtained by FOl and by PKOl 
for the same bursts shows that while many of the estimates agree in two cases 
out of 6 the angles differ by a factor of ~ 3. Furthermore one has to worry 
about the possible angular inhomogeneity discussed earlier. Given these facts 
the narrowness of the E^ distribution is remarkable! 

PKOl have modeled the afterglow emission over a wide range of frequency 
and time for 8 well studied bursts. The model is fitted to the data and the fit 
yields several burst parameters, including the adiabatic energy, Ek, and the jet 
opening angle. The results are also shown on Fig 3. Using their estimate of the 
opening angles 6i PKOl also estimate E^ (see Fig. 3). The PKOl estimates for 
£?-y, differ from the FOl estimates for the same bursts by up to a factor of 8! 

4. The Width of the Energy Distribution 

Both the FOl and PKOl results show that the energy distributions are rather 
narrow. FOl find that the FWHM of the E^ distribution is only a factor of 5. 
PKOl Ek estimates are very narrow, the width is only a factor of 3. The E~^ 
estimates of PKOl are wider, but still narrow with a width of 11. These values 
should be compared with the Eiso distribution that spans a factor of 10"^. This 
motivates us to search for different limits on the width of the energy distribution. 

Following Piran et al (2001a) I turn now to determine the spread of Ek 
using the x-ray afterglow flux. The key factor in this method is the observation 
that the x-ray luminosity at a fixed time after the burst depends almost linearly 
on the kinetic energy of the flow at that time and very weakly if at all on other 
parameters (Kumar 2000; Freedman &; Waxman, 2001). 
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Figure 3. E^^iso (diamonds) from Bloom, et al., (2001); (crosses) from FOl, 
and E-y (circles) and Ek (stars) from PKOl. ^^.(FOl) and E-y{PKQl) for the same 
burst are connected with a vertical dotted line. The Eiso distribution is consistent 
with the average isotropic energy determined by Jimenez et al (2001) for bursts 
with afterglows (dashed-dotted line). Schmidts' (1999) average E^^iao is shown as a 
dashed line. Most E^^iso estimates for burst with afterglow observations are higher 
than Schmidts' 



Consider a fireball with F > 6"^, so that the spherical approximation holds. 

The standard synchrotron fireball model implies that the isotropic equivalent 
luminosity at frequency v above the cooling frequency, at a fixed elapsed time 
since the explosion, is given by (Kumar 2000; Preedman Sz Waxman, 2001): 



Lx = rip 



cLEk 

In 



where dEx/dO, is the kinetic energy per unit solid angle, and 77^ is a constant. 
Lx depends on the energy per unit solid angle in the explosion, and on the 
fractional energy taken up by electrons, eg. p here is the specral index of the 
electron's energy distribution. Since typically w 2 we have an almost linear 
dependence of on Ek- should be measured several hours after the burst 
before any jet break takes place. 

One can attempt to estimate Ek from measurements of L^- However a 
simpler and more robust calculation will be to estimate the width of the energy 
distribTition (JEj^ which in turn is determined by the width of the x-ray luminosity 
distribution. Remarkably, the width of the x-ray luminosity distribution can be 
estimated from the observed x-ray flux distribution. The x-ray afterglow fluxes 
from GRBs have a power law dependence on v and on the observed time t (Piro, 
2000): fi/{t) oc v~l^t~°' with a ~ 1.4 and j3 ~ 0.9. The observed x-ray flux per 
unit frequency, fx^ is related, therefore, to, Lx, the isotropic luminosity of the 
source at redshift, z by: 

Lx{t) = (^^(0(1 + zf-"" = mZ{z) , (3) 

where Z{z) is a weakly varying function of z. For bursts with 0.5 < 2: < 4 and 
with (5 — a K. —0.5 we find az ~ 0.31 (for a cosmology with Q,^ = 0.3 and 
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Figure 4. From Piran et al, (2001) Left: The distribution of X-ray fluxes (2-10 
keV) at t=ll hours after the GRB in 21 afterglows observed by BeppoSAX. The 
arrow marks two bursts for which there is only an upper limit. Right: Likelihood 
contour lines (corresponding to 99%, 90% and 69% confidence levels) in the log(/a;), 
cT/^ plane for the X-ray flux distribution as inferred from 21 GRBs detected by 
BeppoSAX. The maximal likelihood is at log(/a,) = -12.2±0.2 and <t/^ = OASt'll- 



r^A = 0.7). Here and thereafter we denote by ax the standard deviation of the 
log{X), unless noted otherwise. 

Piran et al (2001) use 21 BeppoSAX bursts (Piro, 2000) to determine \og{f^) 
and af^ = 0.43l'Jf for the observed x-ray flux in the 2—10 kev band at 11 
hr after the GRB (For two of the 21 bursts there is only upper limit of 2 x 
lO^'^^ergs/cm^/sec to the x-ray flux. This is consistent with a predicted number 
of 3.5 burst with x-ray afterglows below this limit.). 

If there is no correlation between the microscopic variables eg, e^, p and 
cIEk I dVL and if the x-ray flux does not depend on the redshift we have: 

+ = < <^l. = 4 + 4 = (o-43i:l?)' . (4) 

PKOl and FOl have estimated the jet opening angles obtaining gq^PK^I) 
0.31 ± 0.06 and cre(-FOl) ^ 0.28 ± 0.05 for 8 and 17 bursts respectively. If these 
values are representative for the whole GRB population wc find a marginally 
viable solution within two a errors of < 0.2 (for the PKOl result) and 
(^Ek < 0-27 (for the FOl data); to get a viable solution we had to take both the 
values of gl,^ one a above the mean and the value of ag one a below the mean. 
This result suggests that there is a narrow energy distribution; the FWHM of 
Ek being less than a factor of 5. 

We have argued before that Ek, discussed here, is a rather good estimate 
to Ej-gi the total energy emitted by the "inner engine" . The constancy of Ek 
is another indication for it being a good measure of E^^i- The constancy of Ek 
is also an indication that the assumptions that have lead to Eq. 2 are justified. 
Otherwise it would have been remarkable if starting from different levels of initial 
energy and having different amounts of energy losses the final kinetic energy of 
the afterglow would converge to a constant value. 
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5. Discussion 

There are two striking results: the narrowness of the and Ek distributions 
and the fact that E^ 3Ek (see Fig. 3). There seem to be an intrinsic 
inconsistency between the two results. A narrow Ey distribution with E^ > Ek 
implies (with no fine tuning) that E^ei is narrowly distributed and E^ k, Erei 
(rather than Erei ~ Ek)- The fact that Ek < E^ is also narrowly distributed 
implies that e, the conversion efficiency of relativistic kinetic energy to 7-rays, 
is close to unity and moreover, e itself should be rather narrowly distributed 
(between 70-80%). This last conclusion would be astonishing considering the 
dependence of e on the the distribution of energies and Lorentz factors of the 
different shells. The narrow distribution of also implies that the emitting 
jets are rather homogeneous (otherwise we would expect significant variations 
in E^ even for the same burst when observed from slightly different direction). 
The results are possible but they require a (currently) inexplicable fine tuning. 

There is no simple way out. It is of course possible that we have been misled 
by small number statistics and we have to wait for a larger data set in which we 
expect to find Ek > E^ and a larger spread in E^. 

The current results would have been easier to understand if, for some reason, 
Ek have been underestimated by a factor of 3. In this case we would have 
Erei ~ ~ Ek with e ~ 0.5 and with a reasonable spread in e. Note that at 
least in the PKOl results we already have > (^Ek- 

A second possibility is that in fact Ek > E^ (where E^ is the 7-ray energy 
averaged over the whole GRB jet). However, in our sample angular inhomo- 
geneities resulting in bright spots have increased, according to the patchy shell 
model (Kumar & Piran, 2000) , the observed E-y (in the afterglow sample) above 
the real (average) E^ value. This would results in cte^ > <te^, as observed 
indeed in the PKOl results. However, the same angular inhomogeneities would 
result at times in a decrease in the observed flux. We should observe also bursts 
with E^ < Ek- Such bursts are currently missing in the PKOl sample! Again 
more bursts may resolve this problem. 

In any case these results indicate that one way or another GRB "inner 
engines" are standard candles releasing a rather constant energy. The wide 
distribution of directly and indirectly determined £"-,,^50 results from the wide 
distribution of beaming angles. The fact that GRB engines are "standard" 
engines in terms of their energy output provide a very severe constraint on the 
nature of these enigmatic explosions. For instance, in the collapsar model for 
GRBs the central engine is composed of a black hole (BH) and an accretion disk 
around it (Woosley 1993; Paczynski, 1998; MacFadyen & Woosley, 1999). This 
model has two energy reservoirs which can be tapped to launch a relativistic jet: 
the BH rotation energy and the gravitaional energy of the disk. Our result of 
nearly constant energy in GRBs implies that the mass accretion on to the BH 
plus the possible conversion of rotational energy of the BH to kinetic energy of 
the jet does not vary much from one burst to another in spite of the fact that 
both the disk mass and the BH spin are expected to vary widely in the collapse 
of massive stars. 

I thank P. Kumar, E. Nakar A. Panaitescu and L. Piro for helpful discus- 
sions. This research was partially supported by the US-Israel BSF. 
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